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1.0  INTRODUCTION 


1.1  OBJECTIVE 

In  the  past  several  years,  the  use  of  partially  ionized  gases,  com¬ 
monly  called  plasmas,  in  the  research  and  test  environment  has 
become  commonplace.  These  plasmas  occur  in  such  diverse  applica¬ 
tions  as  arc-heated  wind  tunnels,  lasers,  electron  beam  diagnostics, 
reentry,  and  radiofrequency  (RF)  heaters.  The  manner  of  treatment  of 
the  plasma  state  usually  employed  is  to  express  the  gross  properties  of 
flows  in  terms  familiar  for  un-ionized  gases,  and  at  the  same  time 
accoimt  for  the  observables  peculiar  to  the  plasma  state.  The  bulk  gas 
properties  may  still  be  expressed  as  pressure,  temperature,  density, 
etc. ,  while  the  plasma  properties  are  defined  in  terms  of  electron 
density,  recombination  rate,  etc. 

Because  the  electron-ion  recombination  process  is  necessarily  a 
collisional  one,  it  is  generally  described  through  the  use  of  a  recombi¬ 
nation  rate  coefficient.  The  recombination  rate  coefficient  is  related 
to  the  probability  of  the  occurrence  of  a  successful  recombination  event, 
integrated  over  the  energy  distribution  function  of  the  collision  partners. 
Because  some  cross  sections  (probabilities)  are  very  difficult  and  others 
are  impossible  to  obtain  directly  from  theory,  the  rate  coefficients  are 
generally  obtained  through  experimental  means.  Although  recombination 
rate  coefficient  determinations  have  been  the  subject  of  numerous  inves¬ 
tigations,  the  measurements  are  often  compromised  or  made  more 
difficult  by  the  short  plasma  life  time  of  by  containment  effects. 

The  present  work  is  an  examination  of  the  processes  affecting  the 
electron  number  density  along  the  centerline  of  the  freely  expanding 
plasma.  The  time  rates  of  change  of  the  electron  density  attributable 
to  the  gas  dynamic  and  recombination  mechanisms  are  compared  to 
ascertain  the  utility  of  a  simple  approach  of  determining  the  recombi¬ 
nation  rate  coefficient.  This  approach  requires  only  the  determination 
of  the  electron  number  density  at  points  along  a  streamline  of  the  freely 
expanding  plasma,  as,  for  example,  the  centerline  of  an  arcjet  plume. 

To  be  successful,  the  recombination  rate  must  necessarily  dominate 
the  time  rate  of  change  of  the  electron  number  density  attributable  to 
gas-dynaimic  processes.  The  examination  of  the  relative  magnitude  of 
these  two  rates  constitutes  the  major  portion  of  this  report. 
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1.2  BACKGROUND 

1.2.1  Experimental  Determination  of  Recombination  Rates 

Determining  the  value  of  the  electron-ion  recombination  rate  coeffi¬ 
cient,  a,  for  different  gases  has  been  the  subject  of  numerous  studies 
(Refs.  1  through  11).  In  1896,  Thomson  and  Rutherford  made  the  first 
studies  of  recombination  (Ref.  1)  and  proposed  the  mechanism  of 
recombination  to  explain  the  gradual  decrease  in  the  conductivity  of  a 
gas  following  its  ionization  by  a  burst  of  X-rays.  Various  experimental 
techniques  have  since  been  used  to  determine  the  recombination  rate 
coefficient.  Difficulties  in  experimentally  determining  a  are  due  to  the 
loss  of  electrons  from  a  measurement  volume  by  diffusion  as  well  as 
recombination  (Ref.  1). 

Effective  measurements  of  helium  recombination  rate  coefficients 
Biondi  and  Brown  (Refs.  2  and  3)  utilized  microwaves  to  ionize  the 
sample  gas  as  well  as  a  means  for  measurement.  These  measurements 
depended  upon  obtaining  the  mean  decay  time  of  charge  density  following 
breakdown  in  a  container  of  known  dimensions.  Diffusion  and  wall 
recombination  effects  were  present,  and  it  was  necessary  to  account 
for  them.  Oskam  (Ref.  4)  analyzes  some  of  the  shortcomings  of  this 
approach.  Despite  these  disadvantages,  other  investigators,  including 
Redfield  and  Holt  (Ref.  5),  Oskam  and  Mittelstadlt  (Ref.  6),  and  Chen 
(Refs.  7,  8)  have  used  the  approach  for  ascertaining  the  recombination 
rate  coefficients  for  the  rare  gases  He,  Ne,  Ar,  Kr,  and  Xe.  Other 
investigators  are  cited  in  Ref.  1. 

Another  approach  which  has  been  used  successfully  in  recent  years 
for  both  helium  and  argon  is  that  of  laser  interferometry.  The  laser 
work,  pioneered  by  Gusinow,  Gerardo,  and  Verdeyen  (Ref.  9)  uses  a 
plasma  produced  in  a  laser  cavity.  The  resultant  change  in  the  laser 
frequency  caused  by  the  presence  of  the  electrons  yields  the  necessary 
data  for  the  recombination  rate  coefficient  determination.  Their  proce¬ 
dure  uses  the  refractivity  of  the  plasma  at  optical  wavelengths  and  the 
expression  deduced  by  Gerardo  and  Verdeyen  (Ref,  10).  The  experi¬ 
ment  also  requires  accounting  for  diffusion  effects. 

Another  approach  in  recent  years  has  been  to  use  spectroscopic 
techniques  to  obtain  the  fundamental  data  necessary  to  ascertain  the 
recombination  rate  coefficients.  This  approach,  illustrated  for  example 
by  the  work  of  Shirai  and  Tabei  (Ref.  11),  buHds  upon  an  analytic  tech¬ 
nique  first  suggested  by  Byron,  Stabler,  and  Bortz  (Ref.  12)  and  requires 
detailed  measurement  of  absolute  spectral  intensities  as  well  as  a 
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detailed  calculational  procedure.  Shirai  and  Tabei  used  as  their  plasma 
source  the  freely  expanding  plume  from  a  helium  arcjet,  thus  showing 
the  basic  feasibility  of  the  experimental  approach.  Their  approach 
requires  measurement  of  the  intensities  of  spectral  lines  originating 
from  the  several  excited  states  and  relating  these  to  the  absolute  num¬ 
ber  densities  of  each  energy  level.  The  number  densities  can  then  be 
related  to  the  collisional  cross  sections  and  thence  to  the  recombination 
rate  coefficient. 

An  alternative  possibility,  much  simpler  in  concept,  is  to  determine 
the  free  electron  number  density  at  successive  points  along  the  center- 
line  of  the  plume  from  the  arcjet.  This  number  density,  along  with  an 
accounting  of  the  flow  velocity,  allows  the  average  rate  of  change  of 
the  free  electron  number  density  to  be  calculated  and  thence,  related 
to  the  recombination  rate  coefficient. 

1.2.2  Theoretical  Calculations 

Several  approaches  have  also  been  taken  to  determine  theoretically 
the  recombination  rate  coefficient  for  helium  and  argon.  One  of  the 
methods  is  a  solution  of  the  quasi- equilibrium  equations  for  collision- 
radiative  recombination  as  proposed  by  Bates,  Kingston,  and  McWhirter 
(Ref.  13);  This  method,  descriptive  of  the  physical  processes  involved, 
will  be  discussed  in  more  detail  in  a  subsequent  section  and  was  used  by 
Drawin  (Ref.  14)  to  obtain  recombination  rate  coefficients  for  helium. 
Drawin's  results  are  used  in  the  present  study  for  theoretical  helium 
calculations  coupling  the  gas-dynamic  and  recombination  processes. 

In  a  second  method,  proposed  by  Byron,  Stabler,  and  Bortz  (Ref.  12), 
the  principle  is  used  that  under  equilibrium  conditions  there  exists  a 
pronounced  minimum  in  the  total  rate  of  de -excitation  of  atoms  as  a 
function  of  the  principal  quantum  number  of  the  excited  state.  This  ' 
minimum  limits  the  net  rate  of  three -body  recombination  to  the  rate  of 
de -excitation  of  the  level  at  which  the  minimum  occurs.  Using  this 
simpler  calculational  method  in  the  hydrogenic  approximation,  Wanless 
(Ref.  15),  obtained  argon  recombination  rate  coefficients  in  agreement 
with  the  experimental  values  of  Chen  (Ref.  7).  Since  Wanless'  results 
agree  with  the  more  complex  method  and  cover  a  wide  range  of  electron 
temperatures  and  electron  densities,  they  are  used  in  the  present  study 
for  theoretical  argon  calculations  coupling  the  gas -dynamic  and  recom¬ 
bination  processes. 
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2.0  THEORY 


2.1  INTRODUCTION 

The  description  of  the  free  expansion  of  a  partially  ionized  gas  from 
an  arcjet  involves  two  distinct,  yet  coupled,  physical  phenomena. 
Recombination  processes  are  continually  providing  for  the  adjustment  of 
the  number  densities  of  the  charged  specie  (free  electrons  and  singly 
charged  ions  for  the  purposes  of  this  study).  The  gas-dynamic  processes 
involved  in  the  rapid  expansion  of  the  gas  from  a  reservoir  causes  the 
number  densities  to  decrease  with  increasing  distance  from  the  reser¬ 
voir.  Since  the  description  of  each  phenomenom  involves  quantities 
affected  by  the  other,  a  theoretical  description  of  the  processes  must 
include  considerations  for  superposing  the  two  phenomena.  In  this 
section,  the  equations  describing  recombination  and  gas-dynamic  pro¬ 
cesses  are  discussed,  as  well  as  the  technique  used  for  superposing 
the  two  processes. 


2.2  COLLISIONAL-RADIATIVE  RECOMBINATION 

The  initial  work  on  a  broad  description  of  electron-ion  recombina¬ 
tion  was  done  by  Bates,  Kingston,  and  McWhirter  (Ref.  13)  for  optically 
thin  and  optically  thick  hydrogenic  plasmas.  For  the  case  of  an  opti¬ 
cally  thin  plasma,  they  confined  their  discussion  to  the  recombination 
of  electrons  (denoted  by  e“)  with  bare  nuclei  (denoted  by  N+)  of  charge 
Ze  (where  Z  is  the  level  of  ionization  and  e  is  the  charge  on  an  electron) 
to  form  hydrogen  atoms  or  hydrogenic  ions.  Letting  p,  q  denote  the 
quantum  numbers  of  the  discrete  levels  and  letting  c  represent  the  con¬ 
tinuum,  the  number  densities  of  atoms  or  ions  in  the  indicated  levels 
are  denoted  by  n(p),  n(q).  The  niimber  density  of  free  electrons  are 
designated  by  n(c)  and  the  ions  n(N^). 

The  coUisional  excitation  or  de-excitation  process,  which  can  be 
written 


n(p)  +  e’  -»•  n(q)  e 


is  described  with  K(p,  q)  as  the  rate  coefficient.  Radiative  recombination, 

n(iN’'')  +  e"  n(p)  +  hi/ 


10 


AEDC-TR-76-34 


is  described  by  the  rate  coefficient  /3(p).  Three-body  recombination  into 
some  state  p  is  described  by  the  rate  coefficients  K(c,  p),  and  the  reverse 
ionizing  process  by  the  rate  coefficient  K(p,  c).  The  three-body  reaction 
is  written 

nCN"*)  -I  e*  ^  n(p)  i  e' 

The  spontaneous  radiative  transition  probability  for  the  reaction 

n(p)  ->  n(q)  +  hi-' 

is  denoted  by  A(p,  q). 

With  the  assumptions:  (1)  that  the  electronic  transitions  caused 
by  atom-atom,  atom-ion,  or  ion-ion  collisions  can  be' neglected;  (2) 
boundciry  effects  can  be  neglected;  and  (3)  the  plasma  is  optically  thin 
so  that  all  emitted  radiation  escapes  the  plasma  without  being  absorbed. 
Bates  et  al.  (Ref.  13)  were  able  to  express  the  rate  of  change  of  n(p) 
with  time  by 

=  -n(p)  rn(c)  K(p,q)  -  n(c)  S  K(p.q)  +  v  A(p,q) 

L  ■I’^P  'KP 

+  n(c)  S  n(q)K.(q.p)  +  S  n(q)  A(q,p) 
qi  p  q>  p 

+  n^(c)  [K{c.p)  +  ^(p)] 

The  terms  in  Eq.  (1)  can  be  identified,  respectively,  as  depopu¬ 
lation  caused  by  collisional  ionization,  collisional  excitation,  or  de- 
excitation  from  state  p  to  some  other  q^^  state,  radiative  decay  from 
state  p  to  state  q,  population  of  state  p  caused  by  collisional  excita¬ 
tion  or  de-excitation  from  some  other  q^^  state,  radiative  decay  to 
state  p  from  some  higher  state  q,  three-body  recombination  into 
state  p,  and  radiative  recombination  into  state  p.  Thus,  Eq.  (1)  forms 
the  algebraic  combination  of  all  the  rates  that  add  to  or  subtract  from 
the  number  density  of  atoms  in  level  p.  It  should  be  noted  that,  in  the 
formulation  of  Eq.  (1),  if  a  particular  effect  or  transition  is  not  al¬ 
lowed,  the  corresponding  rate  coefficient  would  be  zero.  If  one  thus 
writes  an  equation  of  the  type  represented  by  Eq.  (1)  for  each  state 
available  to  the  plasma,  one  has  at  hand  the  complete  transient  de¬ 
scription  for  the  plasma  decay.  The  solution  to  that  problem  has  been 
attacked  elsewhere  (Ref.  16). 

However,  following  Bates,  et  al.  (Ref.  1),  assuming  that  a  quasi- 
steady-state  condition  (all  rates  are  negligible  except  dn(l)/dt  and 
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cin(c)/dt)  is  established  almost  instantaneously,  and  utilizing  certain 
eqiiilibria  considerations,  the  system  of  equations  represented  by 
Eq.  (1)  can  be  replaced  by  the  equation 


dn(  1) 
dt 


rn(c)n(N+) 


(2) 


for  single  ionization.  Equation  (2)  expresses  the  time  rate  of  change 
of  the  population  of  the  ground  state,  n(l),  of  the  atom  and  depends 
upon  the  time  rate  of  change  of  all  other  quantum  states  being  negligi¬ 
bly  small  compared  to  that  of  the  ground  state  and  the  free  electron 
density.  In  Eq.  (2),  r  can  be  written  as 


r  =  a 


(3) 


where  a  is  called  the  collisional -radiative  recombination  (CRR)  rate 
coefficient,  and  S  is  the  collisional-radiative  ionization  (CRI)  rate 
coefficient.  Continuing  with  the  results  of  Bates,  et  al. ,  and  substitu¬ 
ting  Eq.  (3)  into  Eq.  (2),  one  obtains 

=  an(c)n(N+)  -  S  n(l)n(N^)  (4) 

d  i 


which  yields  dn(l)/dt  in  terms  of  a  and  S.  For  the  conditions  encoun¬ 
tered  in  this  study,  the  terms  resulting  from  CRI  are  negligibly  small 
and  are  thus  neglected  through  the  remainder  of  this  work.  For  example, 
at  an  electron  temperature  of  10,000®K,  an  atom  density  on  the  order 
of  10^^  cm“3,  and  assuming  one-percent  ionization,  the  CRR  coefficient 
for  helium  is  on  the  order  of  10"!^  cm^/sec  and  the  CRI  coefficient 
for  helium  is  on  the  order  of  10“^®  cm^/sec  (Ref,  14).  This  yields 

rates  on  the  order  of  10^^  cm”^/sec  and  10^^  cm  ^/sec  for  CRR  and 
CRI,  respectively. 

Since  the  plasmas  considered  in  this  work  are  assumed  to  be 
singly  ionized,  the  quasi-steady-state  assumption  yields  the  condition 

dn(l)  dn(c) 

dt  “  “  dt  ^  ' 
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or 


dn(c) 

dt 


-  a  n^(c) 


(6) 


In  order  to  determine  the  recombination  rate  coefficient  in  the  absence 
of  other  effects  and  using  Eq.  (6),  the  electron  number  density  along 
with  the  time  rate  of  change  of  number  density  is  all  that  is  necessary. 


23  GAS  DYNAMICS  AND  RECOMBINATION 


In  the  free  expansion  of  a  singly  ionized  plasma  from  an  arcjet, 
the  change  in  the  electron  number  density  from  one  point  to  another 
along  the  plume  centerline  is  due  not  only  to  recombination,  but  also 
to  the  gas -dynamic  effects.  The  velocity  of  the  flow  can  be  used  in 
relating  time  and  position,  so  that  in  principle,  a  typical  gas-dynamic 
(spatially  dependent)  description  can  be  used  in  determining  the  time 
rate  of  change  of  the  electron  number  density  attributable  to  the  expan¬ 
sion  effects. 


It  has  been  shown  that  if  the  degree  of  ionization  is  small,  conven¬ 
tional  isentropic  gas -dynamic  theory  should  give  a  good  approximation 
for  macroscopic  plasma  properties  (Ref.  17)  in  the  arcjet  plume,  such 
as  pressure,  temperature,  and  density.  Therefore  (Ref.  18), 


My’’-' 


(7) 


where  7  is  the  ratio  of  specific  heats  (7  =  1.67  for  noble  gases),  M  is 
the  Mach  number  (M  =  V/c,  where  V  is  the  velocity  of  the  flow  and  c 
is  the  speed  of  sound),  p  is  the  mass  density  at  some  point  downstream 
of  the  nozzle,  and  Pq  is  the  mass  density  at  stagnation  conditions.  The 
mass  density  is  directly  proportional  to  the  number  of  particles  per 
unit  volume;  therefore,  for  two  points  in  the  flow 

-  ^  M?) 


.  (1  ,  Ml) 
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where  the  numerator  in  the  terms  on  the  left  expresses  the  number 
density  of  heavy  particles  at  each  of  the  two  adjacent  positions.  It  is 
assumed  that  the  forces  acting  on  the  ions  are  principally  gas  dynamic 
in  nature.  Thus  one  can  write 


[  n  +  n{N'*')]  2 
[  n  +  nCN"^)]  j 


1 


y-\ 


(9) 


If  it  is  assumed  that  the  level  of  ionization  is  low, 

n  »  n(N+) 


the  recombination  process  will  not  appreciably  affect  the  atom  number 
density.  Thus,  expanding  Eq.  (9)  algebraically, 

nj  +  n(N+)2  =  «M)ni  +  f(M)n(N+)i 


where 


f(M)  = 


1  + 


-  r 


1  + 


y-  1 


1 


The  terms  involving  the  atom  densities  (the  first  terms  on  each 
side  of  the  equation)  can  be  cancelled  and  one  can  write,  approximately 
for  only  the  gas-dynamic  effects  on  the  ion  density. 


n(N+)2 

n(N+)i 


<10) 


It  is  further  assumed  that  the  ions  are  singly  ionized.  It  can  then  be 
shown  that,  at  the  conditions  of  this  work,  the  Debye  length  is  much 
shorter  than  the  nozzle  exit  diameter  so  that  there  is  no  charge  separa¬ 
tion  and  the  plasma  is  electrically  neutral  locally.  Thus  the  gas-dynamic 
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effects  upon  the  ions  also  describes  the  gas -dynamic  effects  upon  the 
free  electrons.  The  ion  and  electron  number  densities  are  further 
depleted  by  the  effects  of  recombination.  This  can  be  accoimted  for  by 
including  the  integrated  effects  of  the  recombination,  or 


n(c)2 

n(c)j 


1 


M 


M 


1 

n{c) 


; 

1 1, 


2 


aii^(c)dt 


(11) 


K  the  two  points  are  sufficiently  dose  together,  the  integral  may  be 

dosely  approximated  by  the  finite  difference  form  and  one  can  write 

1 


n(c)2 


n(c) 


,  y  —  1  2 


■T  V-1 


n(c) j  a 


•l) 


(12) 


where  the  second  term  is  the  contribution  of  the  recombination  process. 
As  was  described  in  Section  2. 2,  a  is  an  implicit  function  of  electron 
number  density  and  electron  temperature,  and  Eq.  (12)  can  be  solved 
iteratively  for  n(c)2  provided  conditions  at  position  one  and  the  Mg 
are  known. 


With  n(c)2  determined  and  assuming  that  conditions  at  position  one 
are  constant  as  view  from  position  two,  the  time  derivative  of  n(c)2  is 


dn(c)2 

dl 


—  a  n(c)2 


(13) 


Thus  the  rate  of  change  of  n(c)  attributable  to  the  gas -dynamic  effect  is 


1 
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and,  as  before,  the  rate  of  change  attributable  to  recombination  is 


'dn(c)\ 

'  'CRR 


=  -  a  n(c)| 


The  M2  in  Eqs,  (12)  and  (14)  is  required  for  the  solution  of  n(c)2 
and  a  comparison  of  the  rates  attributable  to  the  gas -dynamic  expansion 
process.  A.  B.  White  (Eef.  19)  has  determined  that  the  centerline 
Mach  number  distribution  for  the  arcjet  is  adequately  determined  by 
the  model  proposed  by  Ashkenas  and  Sherman  (Ref,  20).  This  model  is 
used  for  determining  M2  and  dM2/dX  in  this  study  and  is  expressed  by 


M2 


(15) 


v/here  A,  C,  and  X^/D  are  empirical  constants  determined  by  Ashkenas 
and  Sherman  by  comparisons  to  the  method  of  characteristics  solutions. 
The  distance  from  the  exit  plane  of  the  nozzle  is  denoted  by  X2  and  the 
exit  diameter  by  D.  For  a  7  of  1.  67,  Ashkenas  and  Sherman  determined 
A  to  be  3.  26,  C  to  be  0.  3 1,  and  Xq/D  to  be  0.  0075. 

Differentiation  of  Eq.  (15)  with  respect  to  X2  yields  dM2/dX2.  The 
velocity  of  the  flow  at  X2  can  be  written  as 


(16) 


where  c  is  the  speed  of  sound,  and  for  an  isentropic  process  using  the 
equation  of  state. 


dt 


M2(yRTg^)*/2 


(17) 


where  R  is  the  gas  constant,  Tg^  is  the  static  gas  temperature  at  position 
2,  and 


T 


§2 


(18) 


where  Tq  is  the  stagnation  temperature. 
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2.4  METHOD  OF  SOLUTION 


The  analytic  description  of  the  coupled  gas-dynamic  and  recombina¬ 
tion  processes  developed  in  the  preceding  sections  was  solved  to  provide 
am  estimate  of  the  comparison  of  the  electron  number  density  decay  rates 
by  the  two  processes.  The  solution  proceeds  as  a  series  of  initial  value 
problems,  each  solution  providing  the  initial  value  for  the  next  solution. 
Since  the  Ashkenas -Sherman  model  is  valid  only  for  X/D  >  1.0,  this 
value,  1,0,  was  chosen  for  the  starting  point  and  a  time  of  1  x  10“9  sec 
was  arbitrarily  associated  with  this  point. 


To  effect  the  solution,  an  X2  is  chosen  and  from  this,  M2,  Tggi  and 
dX2/dt,  are  calculated  [Eqs.  (15),  (18),  and  (17)].  These  values  are  used 
to  determine  the  average  elapsed  time  from  position  one  to  position  two: 


At  = 


X 


2 


1/2 


Ldi 


(19) 


If  this  average  elapsed  time  is  greater  than  a  predetermined  fraction  of 
the  total  time  elapsed  from  the  start  of  the  calculations,  0. 1  in  this 
work,  X2  is  reduced  and  a  new  average  is  computed.  When  the  average 
time  becomes  less  than  the  aforementioned  fraction,  n(c)2  is  determined, 
Eq.  (12). 

In  this  manner,  for  economy  of  computation,  the  time  increment  is 
kept  small  at  the  beginning  of  the  calculations  when  conditions  are  chang¬ 
ing  rapidly  but  is  allowed  to  increase  as  conditions  change  less  rapidly. 
The  rates  of  decay  of  electron  density  attributable  to  the  gas -dynamic 
effects  and  the  recombination  are  then  calculated  [Eqs.  (14)  and  (16)] 
for  comparison. 


2.5  RECOMBINATION  RATE  COEFFICIENTS 

The  final  parameter  necessary  to  complete  the  calculations  is  the 
recombination  rate  coefficient.  These  data  were  taken  from  Wanless 
(Ref.  15)  for  argon  and  Drawin  (Ref.  14)  for  helium.  Wanless'  resxilts 
for  argon  were  extrapolated  graphically  on  a  log-log  plot  of  the  recombi¬ 
nation  rate  coefficient  versus  electron  temperature  at  constant  electron 
number  density  to  provide  the  rate  coefficients  at  electron  number  densi¬ 
ties  greater  than  10 16  cm~6  and  electron  temperatures  greater  than 
16,000*K.  The  results  of  these  graphical  extrapolations  are  given  in 
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Table  1.  It  was  not  necessary  to  extend  the  helium  recombination  coeffi¬ 
cients  from  brawin  <Ref.  14).  For  use  in  the  computer  calculation,  the 
graphical  data  were  converted  to  tables  and  a  three-point,  two-dimensional 
Lagrangian  interpolation  was  performed  in  the  log-log  plane  to  determine 
the  recombination  coefficient  at  a  specific  value  of  electron  temperature 
and  density. 

Table  1 .  Extrapolated  Values  for  Argon  Recombination  Coefficients 


T  ,  n(c),  a. 


0 

n 

1/cm^ 

3 

cm  /sec 

2,000 

1 

X 

lo'^ 

1.6  X 

10"^ 

2,000 

1 

X 

io'« 

1 .6  X 

10"^ 

4,000 

1 

X 

10^^ 

5.4  X 

lo"^ 

4,000 

1 

X 

10l8 

7.6  X 

10-® 

8,000 

1 

X 

lo'^ 

5.9  X 

O 

1 

l-» 

O 

8;  000 

1 

X 

10^8 

9.4  X 

10-'° 

10,000 

1 

X 

lO^' 

8.9  X 

lo"" 

10,000 

1 

X 

10^8 

1.2  X 

10-'° 

14,000 

1 

X 

10^7 

1.4  X 

10-" 

14,000 

1 

X 

00 

o 

1.6  X 

1 

O 

16,000 

1 

X 

ioi7 

7.8  X 

10-'^ 

16,000 

1 

X 

lo'® 

8.2  X 

10“'^ 

20,000 

1 

X 

10'° 

4.2  X 

10-'3 

20,000 

1 

X 

lO^l 

4.5  X 

10-'® 

20,000 

1 

X 

10^^ 

5.4  X 

10-'® 

20,000 

1 

X 

10^® 

7.4  X 

10-'® 

20,000 

1 

X 

10^^ 

1.4  X 

10-'2 

20,000 

1 

X 

lo'^ 

1.4  X 

10-'2 

20,000 

1 

X 

10^® 

2.2  X 

10-'^ 

20,000 

1 

X 

lo'^ 

3.8  X 

IQ-'^ 

20,000 

1 

X 

10^® 

3.5  X 

io-'2 
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3.0  RESULTS 


3.1  CALCULATED  RESULTS 

The  equations  describing  the  coupled  gas -dynamic  and  recombination 
mechanisms  developed  in  the  preceding  section  were  solved  for  the 
electron  number  density  along  the  centerline  of  a  freely  expanding, 
partially  ionized  helium  and  argon  plume.  The  parameters  necessary 
for  the  calculations  are  the  electron  temperature  profile,  the  initial 
stagnation  gas  and  electron  density  conditions,  and  the  starting  position 
in  the  plume.  The  electron  temperature  profiles  used  are  described  in 
the  sections  presenting  results  of  the  helium  and  argon  calculations, 
and  the  other  initial  conditions,  common  to  both  gases,  are  presented 
in  the  next  section. 

3.1.1  Initial  Conditions 

The  initial  conditions  required  to  obtain  the  solution  to  the  coupled 
equations  are  the  total  gas  temperature,  the  electron  number  density, 
and  the  starting  position  in  the  plume,  X/Dq.  As  indicated  earlier, 
the  Ashkenas -Sherman  model  for  Mach  nimiber  is  valid  only  for  X/D  > 

1.  0,  so  this  value,  1.  0,  was  chosen  for  X/Dq.  Since  the  total  gas  tem¬ 
perature  enters  the  calculations  functionally  as  a  square  root  and  then 
only  in  the  calculation  of  the  gas  velocity,  the  results  of  the  calculations 
will  be  relatively  insensitive  to  the  value  chosen.  Thus  the  same  total 
gas  temperature,  3,  000“K,  was  chosen  for  each  of  the  calculations. 

This  value  is  in  the  range  of  previous  observations  (Refs.  21  through  23) 
for  argon.  The  calculations  were  done  for  a  range  of  electron  number 
densities  and  are  listed  in  Table  2.  Included  in  Table  2  are  the  initial 
values  of  electron  number  density  at  stagnation  conditions,  M  =  0,  and 
at  the  position,  X/D  =  1. 


Table  2.  Initial  Values  of  Electron  Number  Density 


1 

5 

1 

5 


M  =  0  X/D  =  1 


X 

lo" 

-3 

cm 

1.45 

X 

-3 

cm 

X 

lo" 

-3 

cm 

7.24 

X 

-3 

cm 

X 

-3 

cm 

1.45 

X 

-3 

cm 

X 

10» 

-3 

cm 

7.24 

X 

-3 

cm 
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3.1.2  Helium 

The  initial  values  for  the  electron  number  density  listed  in  Table  2 
■were  used  to  determine  the  solution  to  the  coupled  recombination  and 
gas-dynamic  equations  for  helium  for  two  analytic  models  of  the  electron 
temperature  decay.  The  first  model. 


=  3540 


(20) 


is  also  used  later  for  comparative  calculations  for  argon  and  was  chosen 
to  simulate  approximately  the  electron  temperature  characteristics 
observed  in  arcjet  expansions.  The  second  model. 


T 


e 


]i  t 


-0.392 


(21) 


is  similar  in  analytic  form  but  describes  a  much  faster  electron  tempera¬ 
ture  decay.  Equation  (21)  was  prompted  by  observing  that  the  recombi¬ 
nation  rate  coefficient  is  a  decreasing  function  of  electron  temperature. 
Thus,  on  a  comparative  basis,  the  recombination  effects  would  be  much 
larger  with  an  electron  temperature  characteristic  modeled  by  Eq.  (21) 
than  when  modeled  by  Eq.  (20). 


Typical  results  obtained  for  the  helium  expansions  are  shown  in 
Figs.  1,  2,  and  3  for  the  electron  temperature  model  of  Eq.  (20).  The 
left  ordinate  of  each  of  these  figures  is  the  time  rate  of  change  of  the 
free  electron  density  attributable  to  the  gas -dynamic  and  recombination 
effects.  The  abcissa  is  elapsed  time.  The  scale  along  the  top  of  the 
figimes  shows  the  corresponding  X/D  at  the  times  listed  on  the  abcissa. 
Included  on  each  figure,  utilizing  the  right  ordinate,  is  a  plot  of  a  for 
the  conditions  of  that  calculation. 


Figure  1  shows  the  results  for  an  Initial  electron  number  density  of 
1  X  10l7  cm"3,  Fig.  2  the  results  for  an  initial  electron  number  density 
of  5  X  10  cm“3,  and  Fig.  3  the  results  for  an  initial  electron  number 

density  of  1  x  lO^S  cm"^.  As  can  be  seen,  from  these  figures,  recom¬ 
bination  mechanisms  dominate  the  decay  only  for  the  highest  electron 
number  density  presented  and  then  only  early  in  the  flow  development. 
From  the  theory  in  Section  2. 0  it  was  seen  that  the  recombination  rate 
is  quadratically  dependent  upon  the  electron  number  density  Eq.  (6), 
whereas  the  gas -dynamic  rate  is  linearly  dependent  upon  the  number 
density,  Eq.  (14).  Hence,  as  indicated  in  these  figures,  one  would 
expect  the  rapid  rise  in  importance  of  recombination  process  with  increas¬ 
ing  electron  number  density.  When  recombination  h»nomes  dominant. 
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this  quadratic  dependence  causes  the  free-electron  number  density  to 
decay  very  rapidly  until  the  gas -dynamic  process  once  again  dominates 
the  time  rate  of  change  of  the  electron  number  density. 


X/D 

1.00  1.08  1.8  9.6 


10'®  10'^  iC®  10'^ 


t,  sec 

Figure  1.  Time  deveiopment  of  eiectron  number  density  decay  rates 
for  heiium  expansion,  n(c)o  =  1  x  cm-^. 
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The  increasing  rate  of  decay  attributable  to  recombination  during 
the  early  stages  of  the  expansion  (t  <  10 "7  sec)  is  due  to  the  increasing 
recombination  rate  coefficient  with  decreasing  electron  temperature. 
During  this  early  time  the  integrated  effect  of  the  rates  upon  the  electron 
density  is  insignificant  and  it  remains  essentially  constant. 

A  more  succinct  and  informative  way  of  comparing  the  rate  of 
electron  density  decay  attributable  to  the  recombination  and  gas-dynamic 
effects  is  by  plotting  the  ratio  of  the  recombination  rate  to  the  gas -dynamic 
rate  versus  X/D.  In  this  manner,  an  ordinate  greater  than  1.  0  signifies 
a  dominance  of  the  recombination  mechanisms.  These  results  for  the 
electron  temperature  model  of  Eq.  (20)  are  shown  in  Fig.  4.  Thus,  for 
the  conditions  of  these  calculations,  the  recombination  mechanism  is 
dominant  in  the  decay  only  for  the  highest  initial  electron  number  den¬ 
sities  (n(c)Q  >  5  X  10l7  cm“^)  very  close  to  the  nozzle  exit  (X/D  <  1.4). 

A  plot  of  Eq.  (20)  is  also  included  on  Fig.  4  with  the  temperature  scale 
given  as  the  right  ordinate. 

If  one-percent  ionization  is  assumed,  these  conditions  suggest  that 
the  total  atom  number  density  for  such  a  helium  expansion  wo;ild  be 
n^  >  5  X  10^9,  or,  at  3,000°K,  the  total  pressure  would  be  greater  than 
150  torr.  This  is  near  the  upper  range  operating  conditions  for  the 
arcjet  used  for  laboratory  studies  in  this  work.  Figure  5  presents  the 
results  obtained  at  the  same  initial  conditions  but  with  the  electron 
temperature  model  of  Eq.  (21).  Because  the  electron  temperature  for 
these  calculations  falls  much  more  rapidly  with  time  [Eq.  (21)  as 
compared  to  Eq.  (20)]  one  expects  the  recombination  processes  to  have 
a  much  more  dominant  role.  This  is  indeed  the  case  as  is  shown  in 
Fig.  5.  The  recombination  dominates  the  entire  decay  for  each  of  the 
initial  electron  number  densities  in  the  very  early  stages  of  the  decay. 
However,  the  recombination  rate  invariably  becomes  dominated  by  the 
gas-dynamic  rate  a  short  distance  downstream  of  the  nozzle  exit  (X/D  <  2). 

An  interesting  feature  of  these  calculations  is  that  the  ratio  of  the 
recombination  rate  to  gas-dynamic  rate  approaches  the  same  value  as 
a  function  of  X/D  for  a  wide  range  of  initial  number  densities.  This  is 
explained  by  the  increase  in  recombination  rate  with  electron  number 
density  at  the  same  electron  temperature  when  recombination  dominates 
the  gas-dynamic  processes.  Thus  the  electron  number  densities  tend  to 
the  same  X/D  distribution.  If  this  common  distribution  is  reached 
before  the  gas-dynamic  processes  dominate  the  flow,  the  subsequent 
expansion  proceeds  at  the  same  X/D  distribution  of  densities,  and  the 
ratio  of  the  recombination  rate  to  gas -dynamic  rate  is  the  same  for  all 
initial  electron  number  densities. 
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Figure  4.  Comparison  of  electron  number  density  decay  rates  attributable 
to  recombination  and  gas  dynamics,  heiium  expansion,  T. 
modeled  by  slowly  decaying  function. 
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If  the  gas-dynamic  processes  commence  to  dominate  before  the 
common  X/D  distribution  is  reached,  the  subsequent  expansion  proceeds 
at  different  electron  number  densities  and  thus  a  different  ratio  of 
recombination  to  gas -dynamic  rate.  That  this  is  the  case  may  be  seen 
by  comparing  Figs.  4  and  5.  In  Fig.  4,  the  gas -dynamic  rate  dominated 
the  expansion  very  early  in  the  decay.  For  X/D  <  1.4  it  is  seen  that  the 
plots  for  the  two  highest  electron  number  densities  were  convergent. 

For  X/D  >  1.4,  however,  the  spacing  between  the  plotted  lines  tends  to 
be  constant.  In  Fig.  5  it  is  seen  that  the  ratio  of  recombination  rate 
to  gas -dynamic  rate  converged  to  the  same  value  at  X/D  « 1.  4  for  all 
but  the  lowest  electron  density.  Calculations  for  X/D>  1.4  yield  the 
same  X/D  distribution  of  electron  number  density. 

3.1.3  Argon 

The  initial  values  for  the  electron  number  density  listed  in  Table 
2  were  used  to  determine  the  solution  of  the  coupled  recombination  and 
gas -dynamic  equations  for  argon  for  two  models  of  the  electron  tempera¬ 
ture  decay.  The  first  of  these  is  the  same  as  was  used  for  the  helium 
calculations 


=  2  510  t'°-* 


(20) 


and  is  usefvil  for  comparisons  of  the  two  gases.  The  second  model  was 
taken  from  McGregor  and  Brewer  (Ref.  24)  who  presented  experimental 
data  from  an  argon  arcjet.  A  graphical  representation  of  the  data  from 
Ref.  24  is  presented  in  Fig.  6. 

The  results  of  some  of  the  calculations  using  Eq.  (20)  are  shown 
in  Figs.  7,  8,  and  9,  plotted  in  the  same  manner  as  Figs.  1,  2,  and  3. 
A  plot  of  a  versus  time  is  included  in  each  figure. 

Figure  7  shows  the  results  obtained  for  an  initial  electron  number 
density  of  1  x  10 cm"^.  Fig.  8  the  results  for  an  initial  electron 
number  density  of  5  x  10^*^  cm"^,  and  Fig.  9  shows  the  resiRts  for  an 
initial  number  electron  density  of  1  x  10 18  cm“3.  Comparison  of  these 
results  with  those  of  Figs.  1,  2,  and  3  shows  that  the  only  essential 
difference  is  the  X/D  as  a  function  of  time.  This  is  to  be 'expected 
since  the  atomic  weight  of  argon  is  40  compared  to  4  for  helium  and 
causes  the  helium  velocity  to  be  approximately  three  times  that  of 
argon;  hence,  the  helium  travels  about  three  times  as  far  as  argon  in 
the  same  time.  The  lower  argon  velocity  causes  the  gas-dynamic 
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contribution  to  the  time  rate  of  change  of  the  free-electron  number 
density  to  be  smaller  than  for  helium.  There  are  some  differences  in 
the  rate  attributable  to  recombination,  also.  Initially  the  differences 
are  due  to  differences  in  the  recombination  rate  coefficient  for  argon 
and  helium  and,  as  the  flow  develops,  differences  in  the  electron  num¬ 
ber  density  also  contributes  to  the  differences  in  the  recombination 
rate.  As  with  helium,  once  the  gas-dynamic  effects  become  the  domi¬ 
nant  mechanism,  they  remain  so  throughout  the  rest  of  the  decay. 


X/D 

Figure  6.  Centerline  electron  temperature  profile,  argon  expansion 
(from  Ref.  24). 
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Figure  8.  Time  development  of  electron  number  density  decay  rates 
for  argon  expansion,  n{c)o  =  5  x  10”  cm-3. 
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Figure  9.  Time  development  of  electron  number  density  decay  rates 
for  argon  expansion,  n(c)o  =  1  x  10“'®  cm  ®. 
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The  results  of  these  calculations  as  well  as  the  other  initial  condi¬ 
tions  are  summarized  in  Fig.  10.  The  ordinate  is  the  ratio  of  the 
recombination  rate  of  change  of  the  electron  number  density  to  the  . 
gas -dynamic  rate  of  change,  and  the  abcissa  is  the  axial  distance  from 
the  nozzle  exit,  X/D.  Comparisons  of  these  results  to  those  shown  in 
Fig.  4  support  the  previous  observations,  with  the  adjustments  appro¬ 
priate  to  the  different  atomic  weight.  As  with  helium,  the  gas-dynamic 
processes  dominate  the  expansion  after  the  expansion  reaches  X/D  «  1.  6, 

The  results  obtained  from  calculations  based  upon  the  electron  tem¬ 
perature  model  reported  in  the  literature  (Eef.  24)  and  shown  in  Pig.  6, 
are  shown  in  Fig.  11.  For  this  profile  also,  as  with  the  analytic  model, 
Eq.  (20),  the  gas -dynamic  mechanisms  dominate  the  decay  for  X/D  >1.6, 

The  calculations  for  helium  and  argon  suggest  strongly  that,  because 
of  the  contribution  of  the  gas -dynamic  processes,  direct  observation  of 
the  free-electron  density  is  not  sufficient  for  recombination  rate  coeffi¬ 
cient  determinations  in  the  expanding  arcjet  plume.  For  the  conditions 
used  in  this  study,  which  are  typical  of  those  encountered  in  the 
laboratory,  observation  of  the  restilts  shows  that  whenever  the  recombi¬ 
nation  mechanisms  do  dominate  the  decay,  their  nature  is  such  that  the 
electron  density  decays  rapidly  until  the  gas-dynamic  processes 
become  dominant. 


3.2  EXPERIMENTAL  DATA 

The  arcjet  has  been  used  in  continuing  research  work  at  AEDC  to 
investigate  various  excitation  mechanisms.  In  the  course  of  this  work, 
heretofore  unpublished  data  have  been  taken  which  can  be  used  for 
comparisons  with  the  results  of  the  previous  sections.  In  this  section 
the  procedures  for  obtaining  the  experimental  data  and  the  analysis  wUl 
be  discussed. 

3.2.1  Apparatus  and  Procedure 

The  plasma  jet  used  is  of  the  Gerdien  type  and  one  model  has  been 
described  in  detail  in  Refs.  25  and  26.  In  the  present  study  the  plasma 
jet  and  other  apparatus  are  most  similar  to  that  described  in  Ref.  26. 
Briefly,  the  arcjet  consists  of  a  hollow  body  into  which  the  working  gas 
is  introduced.  The  plasma  is  produced  in  the  arc  region  between  a 
coaxial,  water-cooled,  tungsten  cathode  and  a  water-cooled,  constant- 
area,  copper  nozzle  anode.  The  plasma  exhausts  through  this  nozzle 
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Figure  10.  Comparison  of  electron  number  density  decay  rates  attributable 
to  recombination  and  gas  dynamics,  argon  expansion,  T, 
modeled  by  slowly  decaying  function. 
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Figure  11.  Comparison  of  electron  number  density  decay  rates  attributable 
to  recombination  and  gas  dynamics,  argon  expansion,  T, 
modeled  from  Ref.  24. 
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into  the  test  cell.  The  pressure  in  the  test  cell  is  maintained  by  two 
large  mechanical  vacuum  pumps.  A  quartz  window  is  mounted  in  the 
test  cell  for  spectroscopic  observations.  In  this  study  the  arcjet 
nozzle  had  a  diameter  of  0.  635  cm.  A  schematic  representation  of  the 
test  apparatus  (not  to  scale),  is  shown  in  Fig.  12. 


The  chamber  and  lines  were  evacuated  before  arcjet  operation  in 
order  to  ensure  minimum  contamination.  The  argon  gas  supply  was  of 
99.  9-percent  purity  (with  oxygen  1  ppm,  nitrogen  1  ppm,  hydrogen 
1  ppm,  and  moisture  1  ppm)  was  supplied  to  the  gas  lines  through  a 
regulator  adjusted  to  maintain  a  line  pressure  of  5,  750  torr.  The 
power  supply  to  the  arcjet  was  a  motor  generator  providing  an  open  line 
voltage  of  85  v.  The  argon  gas  was  introduced  into  the  arcjet  chamber 
and  the  plasma  jet  was  started  with  an  RF  discharge  of  short  duration 
which  is  used  to  break  down  the  gas  and  then  is  shut  off  after  the  arc  is 
started.  Stable  operating  conditions  were  as  follows: 


Chamber  pressure 
Gas  flow 

Operating  amperage 
Operating  voltage 
Cell  pressure 


142.  2  torr 
3.  0  gm/sec 
200  amp 
25  V 
0.  4  torr 


A  1-m  Jerrell-Ash  spectrometer  was  placed  so  that  the  entrance 
slit  (slit  width  of  18  ju )  was  aligned  with  the  centerline  of  the  nozzle 
in  order  to  measure  the  intensities  of  spectral  lines  throughout  the 
radial  extent  of  the  plasma,  the  spectrometer  field  of  view  was  scanned 
across  the  plasma.  This  was  accomplished  by  an  optical  system  con¬ 
sisting  of  a  converging  lens  with  a  focal  length  of  12  in.  and  a  plane 
mirror.  The  mirror  was  rotated  at  an  angular  velocity  of  72  deg/min 
to  provide  the  scan  of  the  field  of  view.  The  lens  focused  the  light 
from  the  plasma  onto  the  slit  of  the  spectrometer.  The  detector  was  an 
RCA  1P28  photomultiplier  tube  which  was  connected  to  a  PhUbrick 
(model  P2)  amplifier.  The  gain  on  the  amplifier  was  set  at  100.  The 
data  were  recorded  on  an  x-y  recorder. 


After  each  scan  the  arcjet  was  lowered  0,  635  cm  and  another  scan 
was  made.  This  procedure  was  continued  until  measurements  of  the 
intensity  (radiated  power  per  unit  area  per  unit  solid  angle)  of  the 
selected  spectral  lines  had  been  made  at  0.  635-cm  intervals  along  the 
plume  axis  over  arrange  of  travel  3  <  X/D  <  6.  The  4 158- A  spectral 
line  and  the  6032-A  spectral  line  were  selected  for  the  radial  scans 
because  reliable  transition  probabilities  are  available  (Ref.  27)  for 
these  lines. 
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A  calibrated  tungsten  filament  lamp  was  used  for  calibration  of 
the  spectrometer  and  optical  system  as  shown  in  Fig.  12.  Table  3  lists 
the  calibration  scale  factors  for  the  two  wavelengths  (4158  A  and  6032  A) 
that  were  used  in  this  study.  With  these  calibration  values  (which  take 
into  account  the  difference  in  using  the  tungsten  source  as  compared  to 
the  plasma  source)  the  numbers  read  from  the  traces  on  the  x-y 
recorder  can  be  given  physically  meaningful  units.  Thus  the  raw  data 
are  related  to  a  true  intensity  of  w/ster-cm^. 


Figure  12.  Research  cell,  plasma  jet,  and  optics  installation. 


Table  3.  Intensity  Calibration  Factors 


X/D 

Intensity, 

2 

w/ster-cm  /div 

4158 

O 

A 

6032 

0 

A 

3.0 

7.12  X 

10-‘ 

3.63  X 

10-6 

4.0 

7.12  X 

10-® 

3.73  X 

10-6 

5.0 

2.85  X 

10~^ 

1.87  X 

10-6 

36 


AEDC-TR-76-34 


3.2.2  Data  Analysis 

Figures  13  through  18  are  traces  from  the  x-y  recorder  of  the 
radiation  intensities  for  both  the  4158-A  line  and  the  6032-A  line  that 
were  used  in  this  study  for  X/D  =  3,  4,  and  5.  The  plume  was 
scanned  at  least  twice  for  each  spectral  line  at  each  axial  position  to 
verify  arcjet  stability  and  symmetry.  The  center  of  each  peak 
corresponds  to  the  location  of  the  centerline  of  the  plume,  and  the 
scans  were  examined  for  repeatiability  and  symmetry.  In  each  case 
symmetry  was  satisfactory  and  the  peak  height  was  repeatable  within 
5  percent.  The  variations  in  the  two  scans  were  thus  judged  to  be 
insignificant,  and  radiation  intensity  data  as  a  function  of  distance 
from  the  centerline  were  obtained  from  one  side  of  one  scan  only. 

These  data  were  then  converted  to  intensity  units  by  the  appropriate 
calibration  factor  given  in  Table  3. 

The  spectral  line  intensity  measurements  were  used  to  obtain  the 
local  emission  coefficients  (radiated  power  per  unit  solid  angle)  by 
use  of  a  computer  code  designed  to  effect  the  Abel  inversion  integral. 
The  approach  utilized  in  the  code  was,  via  a  least-squares  spline  fit, 
to  model  the  original  intensity  profiles  with  a  series  of  polynomials 
of  even  powers  of  distance.  The  even  function  is  chosen  to  ensure 
cylindrical  symmetry,  a  requirement  of  the  Abel  inversion  for  this 
application.  With  the  intensity  modeled  by  the  appropriate  series  of 
polynomials,  the  inversion  is  analytic  and  immediate.  Figure  19 
presents  representative  results  of  these  inversions  with  the  emission 
coefficient  normalized  by  the  centerline  value  for  data  at  one  of  the 
points  of  measurement  in  this  study. 

Once  the  emission  coefficients  are  calculated  the  electron  tempera¬ 
ture  is  determined  with  the  assumption  of  local  thermodynamic 
equilibrium  (Ref.  28).  The  emission  coefficient  is  divided  by  the 
statistical  weight  of  the  emitting  state,  the  transition  probability  for 
the  transition,  and  the  frequency  of  the  emitted  radiation.  The 
logarithm  of  this  value  is  then  plotted  against  the  energy  of  the  radiat¬ 
ing  state.  These  plots  give  nearly  straight  lines  of  slope  -0.4343/kTe, 
yielding  immediately  the  electron  temperature.  In  Fig.  20  is  a  plot 
of  the  electron  temperature  as  determined  from  the  emission 
coefficients  along  the  centerline.  The  bars  indicate  the  observed 
variation  in  the  calculated  electron  temperature  from  variations  in  the 
fitted  intensity  data  within  the  previously  mentioned  5-percent  variation. 
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Spatial  Displacement  across  the  Plume 


Figure  13.  Radial  plasma  jet  radiation  intensity  scans,  X/D  =  3, 
X  =  4158  A. 


Spatial  Displacement  across  the  Plume 


Figure  14.  Radial  plasma  jet  radiation  intensity  scans,  X/D  =  4, 
X  =  4158  A. 
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Spatial  Displacement  across  the  Plume 


Figure  15.  Radial  plasma  jet  radiation  intensity  scans,  X/D  =  5, 
X  =  4158  A. 


Figure  16.  Radial  plasma  jet  radiation  intensity  scans,  X/D  =  3, 
X  =  6032  A. 
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Spatial  Displacement  across  the  Plume 


Figure  17.  Radial  plasma  jet  radiation  intensity  scans,  X/D  =  4, 
X  =  6032  A. 


Figure  18.  Radial  plasma  jet  radiation  intensity  scans,  X/D  =  5, 
X  =  6032  A. 
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Since  the  gas  is  assumed  to  be  singly  ionized,  the  ion  number 
density  and  electron  number  density  are  equal.  Using  the  above  deter¬ 
mined  Te  the  n(c)  can  be  determined  from  the  principle  of  mass  action 
(Saha’s  equation)  by  determining  the  n(c)  in  equilibrium  with  an  excited 
state  at  the  measured  electron  temperature  (Ref.  28).  These  resialts 
are  plotted  in  Fig.  21  for  the  electron  number  density  along  the 
centerline  versus  X/D.  The  bars  in  Fig.  21  have  the  same  meaning  as 
in  Fig.  20.  The  three  data  points  indicated  in  Fig.  21  can  be  connected 
with  a  straight  line  on  the  log -log  plot.  If  gas -dynamic  effects  are  the 
dominant  feature  of  the  electron  density  decay,  potential  flow  theory 
demands  that  the  n(c)  decays  according  to  an  inverse  square  law,  i.  e.  , 
on  a  log-log  plot  such  as  Fig.  21,  a  straight  line  of  slope  -2  will 
result.  If  there  are  other  completing  phenomena  for  the  electron  den¬ 
sity  decay,  such  as  the  recombination  process,  these  will  be  in  addi¬ 
tion  to  the  gas -dynamic  effects,  and  the  log -log  plot  will  necessarily 
have  a  slope  steeper  than  -2.  The  straight  line  slopes  obtained  from 
Fig,  21  have  a  median  value  of  -1.  8  with  lower  and  upper  bounds  of 
-1.45  and  -2.  15,  respectively.  Thus,  one  is  led  to  the  conclusion 
that  the  gas -dynamic  effect  is  the  principal  contributor  to  the  electron 
number  density  decay. 

An  alternate  way  to  examine  the  relative  effect  of  recombination 
and  gas -dynamic  expansion  on  the  local  electron  number  density  is  as 
follows.  If  the  physical  effects  of  the  recombination  process  are  to  be 
observed  in  the  experimental  data,  there  should  be  a  noticeable  change 
in  the  ratio  of  the  electron  number  density  to  atom-number  density  as 
the  distance  from  the  nozzle  exit  increases.  The  ratio  L  given  by 

L  =  n{c)/n 


is  given  in  Table  4  for  each  of  the  three  X/D  positions,  3,  4,  and  5. 

For  these  determinations  the  total  number  density  was  obtained  from 
potential  flow  theory  using  stagnation  conditions  determined  by 
applying  Rayleigh  heating  relationships  to  the  arcjet  constant  area  anode 
(Ref.  23).  Table  4  shows  that,  for  each  of  the  three  positions,  L  is 
relatively  constant  with  a  variation  of  only  approximately  8  percent. 

The  almost  constant  value  of  L  implies  that  the  degree  to  which 
recombination  is  affecting  the  electron  density  must  be  very  small 
compared  to  the  gas  expansion. 
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Table  4.  Degree  of  Ionization  at  Appropriate  X/D 

n,  n(c), 


X/D 

—  j 

cm 

—  J 

cm 

L 

3 

2.42  x  10^^ 

2.39  X  10^^ 

9.88  X 

10"^ 

4 

1.59  X  10^^ 

1.45  X  10^^ 

9.12  X 

lO"^ 

5 

1.02  X  10^^ 

9.47  X  10^^ 

9.29  X 

10"^ 

3.2.3  Comparisons  with  Theory 

The  results  presented  above  have  been  consistent  with  the  theore¬ 
tical  predictions  and  analysis  presented  in  Section  3.  1.  However,  the 
experimental  data  presented  here  were  at  a  higher  initial  total  gas 
temperature  (4,  890®K),  whereas  the  electron  number  density  (3.45  x 
1017  cm“^)  was  within  the  range  of  the  conditions  listed  in  Section  3. 1. 
The  higher  gas  temperature  will  affect  principally  the  Mach  number 
of  the  expansion  and  should  not  affect  the  general  conclusion  of 
Section  3.  1  that  gas -dynamic  effects  dominate  the  centerline  rate  of 
decay  of  n(c)  for  X/D  >  2.  However  for  completeness,  the  coupled 
recombination  and  gas -dynamic  program  was  used  with  the  above 
initial  conditions  in  order  to  determine  if  the  resijlts  of  theory  and 
experiment  were  indeed  consistent.  In  Fig.  22  is  a  plot  of  the  time 
rate  of  change  of  the  free-electron  number  density  attributable  to 
gas-dynamic  effects  and  recombination  effects  for  the  initial  conditions 
of  the  experimental  data  versus  time,  and  using  the  computation  method 
described  in  Section  3.  1. 

The  coUisional -radiative  recombination  rate  does  dominate  in  the 
region  X/D  <  2.  However,  in  the  practical  sense  this  portion  is 
inaccessible  because  of  the  blown  arc  which  protrudes  to  about  X/D  =  2 
and  causes  the  assumptions  of  a  field  free  expansion  to  be  invalidated. 

In  the  region  (3.  0  <  X/D  <  6.  0)  where  analysis  of  the  data  is  justifiable 
the  gas-dynamic  rate  dominates,  which  is  consistent  with  the  theoretical 
work  of  Section  2.  1  and  the  analysis  of  the  experimental  data  presented 
in  this  Section. 


4.0  SUMMARY  AND  CONCLUDING  REMARKS 


This  work  has  provided  a  study  of  the  electron  number  densities 
and  recombination  rates  to  be  expected  along  the  centerline  of  the  freely 
expanding  plume  from  an  argon  or  helium  arcjet.  The  purpose  of  the 
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study  was  to  examine  the  relative  rates  of  decay  of  the  free-electron 
number  density  attributable  to  the  recombination  processes  and  the 
gas-dynamic  processes.  For  this  purpose  the  equations  describing  the 
recombination  process  were  superposed  with  the  equations  describing 
the  free  isentropic  expansion  of  an  ideal  gas.  In  this  manner,  calcvila- 
tions  comparing  the  rate  of  change  of  the  free-electron  number  density 
caused  by  recombination  to  the  rate  of  change  caused  by  the  gas -dynamic 
expansion  were  made  for  a  range  of  conditions  similar  to  those  antici¬ 
pated  for  a  typical  arcjet  plume. 

In  the  range  of  stagnation  conditions  studied  (3,  000®K  <  total  gas 
temperature  <  4,  GOO^Kj  1  x  10^"^  cm~3  <  electron  number  density  < 

5.  00  X  10^®  cm"3;  and  1,  200®K  <  electron  temperature  <  20,  000®K)  the 
gas-dynamic  rate  of  decay  of  the  free-electron  number  density  generally 
dominated  the  recombination  rate.  It  was  further  found  that  those 
conditions  for  which  the  recombination  rate  did  dominate  over  the 
gas-dynamic  rate  of  decay  of  the  free-electron  number  density  required 
electron  densities  which  are  in  the  range  of,  but  are  generally  greater 
than,  those  normally  achieved  in  typical  arcjet  experiments.  Further, 
this  domination  by  recombination  processes  occurs  only  for  a  short 
distance  downstream  of  the  nozzle  exit,  X/D  <  2.  Analysis  of  available 
experimental  data  confirmed  these  calcxilated  results.  This  data 
analysis  showed  that  the  electron  number  density  decay  along  the 
centerline  of  the  plume  was  functionally  an  inverse  square  dependence, 
as  is  required  by  potential  flow  theory,  and  further,  that  the  degree  of 
ionization  was  very  nearly  constant  as  a  function  of  axial  distance 
from  the  arcjet  nozzle.  This  latter  result  is  fiirther  evidence  that  the 
recombination  process  did  not  affect  the  electron  density  appreciably 
in  the  region  studied,  as  it  would  have  if  it  were  an  important  feature 
of  the  decay. 

It  is  thus  concluded  that  determinations  of  the  electron  number 
density  in  the  expanding  pliune  of  the  arcjet  will  reflect  the  effect  of 
the  gas -dynamic  expansion  rather  than  the  contribution  of  the  recombi¬ 
nation  mechanisms.  Hence,  the  arcjet  plume  cannot  be  used  for  direct 
determination  of  the  recombination  rate  coefficient  from  the  electron 
number  density.  Rather,  a  previous  approach,  described  in  Ref.  11, 
which  requires  detailed  spectral  intensity  and/or  absorption  measure¬ 
ments  to  determine  the  excited  state  distribution  function  should  be 
used.  These  detailed  intensity  measurements  will  not  only  provide 
means  to  determine  the  recombination  rate  coefficients,  but  also  the 
means  to  examine  the  fully  transient  plasma  decay  problem  (Ref.  16). 
These  studies  could  be  accomplished  without  recourse  to  the  simplfying 
quasi-steady-state  assumption  used  here. 
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NOMENCLATURE 


Empirical  constant 

Transition  probability  from  p  to  q  sec“^ 

Empirical  constant 

Speed  of  sound,  cm/sec 

Diameter  of  nozzle,  cm 

Charge  on  electron,  4.  80286  x  10“^^  stat 
coulomb 

Electron 

Flanck's  constant,  6.  62517  x  lO'^*^  ergs-sec 

Rate  coefficient  for  three-body  recombination 
into  quantum  state  p,  cm3 /sec 

Rate  coefficient  for  coUisional  ionization  by 
electrons  from  quantum  state  p,  cm3 /sec 

Rate  coefficient  for  excitation  or  de-excitation 
from  level  p  to  level  q,  cm3 /sec 

Boltzmann's  constant,  8.6164  x  10“^  eV/®K 

Level  of  ionization 

Mach  number 

Singly  ionized  atom 

Neutral  atom  number  density,  cm"^ 

Electron  number  density,  cm" 3 
Ion  number  density,  cm"^ 

Number  density  of  atoms  in  quantum  state  p,  cm"3 

Symbolic  for  set  of  quantum  numbers  defining 
energy  levels 

Gas  constant,  erg/(g-mole)/“K 

Collisional-radiative  ionization  coefficient, 
cm^/sec 

Electron  temperature,  “K 
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Tg 
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SUBSCRIPTS 
1.  2 

CRR 

G 

o 


Gas  temperature,  °K 
Time,  sec 

Average  elapsed  time,  sec 

Velocity  of  flow,  cm/sec 

Distance  from  nozzle  along  centerline,  cm 

Nondimensional  geometric  constant 

Refers  to  the  level  of  ionization 

Collisional-radiative  recombination  rate 
coefficient,  cm^/sec 

Rate  coefficient  for  radiative  recombination 
into  quantum  state  p,  cm^/sec 

Collisional-radiative  decay  coefficient,  cm^/sec 

Ratio  of  specific  heats 

Wavelength,  cm 

Frequency,  1/sec 

n 

Mass  density,  g/cm^ 

Refers  to  two  points  in  the  flow  along  the 
centerline 

Refers  to  collisional-radiative  recombination 
contributions  only 

Refers  to  gas -dynamic  contributions  only 
Refers  to,  stagnation  values  or  starting  position 


52 


